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1.  Introduction 


Scintillation  caused  by  density  fluctuations  (irregularities)  in  the  ionosphere 
results  in  outages  of  the  communication  and  navigation  systems  that  depend  on  trans- 
ionospheric  radio  links.  The  outages  can  be  particularly  severe  at  low  magnetic  latitudes 
in  the  post-sunset  hours,  when  scintillation  is  caused  by  the  occurrence  of  so-called 
plasma  ‘bubbles’  (density  depleted  regions)  that  are  formed  as  a  result  of  the  generalized 
Rayleigh-Taylor  instability.  Requirements  for  forecasting  ionospheric  scintillation  are 
stated  in  the  Air  Force  Space  Command  Aerospace  Weather  Information  Requirements 
document,  dated  02  April  2002.  The  required  lead-time  for  scintillation  forecasting 
varies  from  1  hour  to  beyond  120  hours  depending  upon  the  individual  unit  and  planning 
cycle  supported.  Current  operational  capabilities  for  scintillation  forecasting  are  limited 
to  climatological  assessments  from  WideBand  ionospheric  scintillation  MODel 
(WBMOD)  and  persistence  techniques  (eastward  drift  and  decay  of  observed  scintillation 
structures)  as  used  in  the  Operational  Space  Environment  Network  Display  (OpSEND) 
Ultra  High  Frequency  Scintillation  Maps.  These  state-of-the-art  operational  maps  provide 
only  a  current  specification  and  a  limited,  two-hour,  forecast  of  ionospheric  scintillation. 

The  objective  of  the  present  task  is  to  develop  a  scheme  for  forecasting 
scintillation  in  the  post-sunset  equatorial  ionosphere  using  NPOESS  or  NPOESS-like 
data.  The  plan  is  to  accomplish  the  task  in  two  phases.  This  report  deals  with  the  recently 
completed  phase  1  of  the  task. 

The  peak  value  of  the  E  x  B  vertical  plasma  drift  velocity  due  to  the  post-sunset, 
enhanced  eastward  electric  field  E  is  now  recognized  as  the  most  important  parameter  for 
determining  both  the  intensity  and  the  latitudinal  distribution  of  equatorial  scintillation 
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{Fejer  et  al.,  1999;  Whalen,  2001].  Estimation  of  the  peak  value  of  the  E  x  B  vertical 
plasma  drift  velocity  is  therefore  critical  to  any  scintillation-forecasting  scheme.  In  our 
proposed  scheme,  the  peak  value  of  the  vertical  plasma  drift  velocity  will  be  estimated 
from  measurements  of  the  Appleton  anomaly  parameters,  such  as  the  magnitude  of  the 
anomaly  crest  and  its  latitudinal  width,  peak-to-valley  ratio  of  the  anomaly  Total  Electron 
Content  (TEC),  etc.  It  will  then  be  used  to  drive  the  Air  Force  Research  Laboratory 
(AFRL)  model  [Retterer,  1999;  Retterer  et  al,  2002]  for  simulation  of  plasma  ‘bubble’ 
formation  and  calculation  of  amplitude  scintillation  index  (S4),  leveraging  the 
developmental  work  within  the  C/NOFS  program.  The  scintillation-forecasting  scheme 
will  thus  be  based  on  the  knowledge  of  the  parameters  characterizing  the  anomaly.  These 
parameters  can  be  derived  from  the  ultraviolet  (UV)  imagery  data  of  the  anomaly  region 
acquired  by  NPOESS  sensors  and  NPOESS-heritage  sensors;  e.g.  GUVI  on  TIMED  and 
SSUSI  on  DMSP.  Hence,  NPOESS  and  NPOESS-like  data  will,  in  principle,  be  capable 
of  being  used  to  forecast  the  occurrence  of  scintillation  and  the  magnitude  of  the 
amplitude  scintillation  index  S4. 

As  a  first  step  towards  developing  the  scintillation-forecasting  scheme,  a 
theoretical  study  was  undertaken  to  derive  and  validate  possible  quantitative  relationships 
between  the  peak  value  of  the  vertical  plasma  drift  velocity  and  one  or  more  of  the 
various  Appleton  anomaly  parameters.  The  relationship  can  then  be  used  to  estimate  the 
peak  value  of  the  vertical  plasma  drift  velocity  from  measurements  of  the  anomaly 
parameter.  The  AFRL  low-latitude  ionosphere  model  [Retterer,  1999;  Retterer  et  al, 
2002]  was  used  for  the  study.  The  scientific  rationale  for  the  feasibility  of  such  a 
scintillation-forecasting  scheme  is  discussed  in  more  detail  in  the  following  section. 
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2.  Scientiflc  Rationale 


An  interesting  feature  of  the  ionosphere  at  equatorial  latitudes,  which  makes  it 
quite  different  from  that  at  other  latitudes,  is  the  so-called  Appleton  anomaly.  During 
daytime,  the  eastward  E-region  dynamo  electric  field  (E)  from  off-equatorial  regions 
maps  along  the  Earth’s  magnetic  field  (B)  into  F-region  heights  above  the  magnetic 
equator.  The  resulting  E  x  B  drift  causes  the  F-region  plasma  to  move  vertically  upward 
at  the  magnetic  equator.  The  uplifted  plasma  then  diffuses  along  B  due  to  gravity  and 
pressure  gradient  forces.  This  results  in  the  formation  of  the  Appleton  anomaly  with 
minimum  F-region  ionization  density  at  the  magnetic  equator  and  maximum  at  the  two 
crests  at  about  ±14  to  ±20  degrees  magnetic  latitude  depending  upon  the  peak  value  of 
the  vertical  plasma  velocity.  Near  sunset,  the  E-region  density  and  the  E-region  dynamo 
electric  field  decrease  and  the  Appleton  anomaly  starts  to  diminish.  However,  after  sunset 
when  the  E  layer  of  the  equatorial  plasma  virtually  disappears  due  to  rapid  recombination 
of  ions  and  electrons,  the  F-region  dynamo,  in  the  presence  of  the  sunset  induced 
conductivity  gradient,  regenerates  the  eastward  electric  field  with  significantly  enhanced 
value.  The  vertical  uplift  of  the  plasma  due  to  the  F-region  dynamo  electric  field  and  the 
subsequent  diffusion  process,  as  described  above,  result  in  the  resurgence  of  the  Appleton 
anomaly  in  the  post-sunset  period.  At  the  same  time,  since  the  recombination  rates 
decrease  with  increasing  altitude,  a  sharp  density  gradient  develops  on  the  bottom  side  of 
the  F  layer,  and  the  plasma  on  the  bottom  side  of  the  F  layer  becomes  unstable  due  to  the 
combined  effects  of  the  gravity,  the  eastward  electric  field  and  the  vertically  downward 
neutral  wind  velocity  (if  there  is  any).  The  relevant  plasma  instability  is  known  as  the 
generalized  Rayleigh-Taylor  (GRT)  instability.  Subsequent  nonlinear  development  of  the 
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instability  leads  to  the  occurrence  of  ‘bubbles’,  which  cause  scintillation.  In  the  most 
elementary,  yet  quite  informative,  description  of  the  GRT  instability,  the  linear  growth 
rate  (y)  is  given  by  [see  e.g.,  Kelley,  1989;  Basu,  1997,  1998,  2002] 

Y  =  (g/Vin  +  VE  +  VN)/Ln  (1) 

Here  g  is  the  acceleration  due  to  gravity,  Vjn  is  the  ion-neutral  collision  frequency,  Ve  is 
the  peak  value  of  the  vertically  upward  E  x  B  velocity  of  the  plasma,  Vn  is  the  vertically 
downward  neutral  wind  velocity  and  Ln  is  scale  length  of  the  density  gradient  in  the 
vertically  upward  direction.  As  is  evident  from  Eq.  (1),  the  peak  value  of  the  vertical 
plasma  velocity  Ve  directly  contributes  to  the  linear  growth  rate.  It  also  enhances  the 
gravitational  term  (g/vin)  in  Eq.  (1)  by  moving  the  plasma  upward  to  higher  altitudes 
where  Vin  is  smaller.  Consequently,  it  plays  an  important  role  in  the  linear  growth  of  the 
GRT  instability  and,  hence,  in  the  subsequent  formation  of  ‘bubbles’  by  nonlinear 
processes.  Furthermore,  scintillation  is  most  disruptive  [Whalen,  2000]  when  a  ‘bubble’ 
encounters  the  maximum  electron  density  of  the  Appleton  anomaly,  which,  as  described 
earlier,  is  also  formed  due  to  the  vertical  drift  of  the  plasma,  albeit  by  a  different  physical 
mechanism.  The  post-sunset  vertical  drift  of  the  plasma  does  not  attain  sufficient 
magnitude  on  some  days  to  cause  the  evening  equatorial  anomaly  and  the  scintillation. 

The  preceding  discussion  attests  to  the  fact  that  the  peak  value  of  the  post-sunset 
vertical  plasma  velocity  plays  an  important  role  both  in  the  formation  of  the  Appleton 
anomaly  and  in  the  occurrence  of  scintillation.  Previous  theoretical  calculations  of  the 
ambient  plasma  density  distribution  in  the  nighttime  equatorial  ionosphere  have  shown 
that  the  structure  of  the  Appleton  anomaly  indeed  depends  on  the  peak  value  of  the  post¬ 
sunset  plasma  drift  velocity  [Anderson,  1973].  In  addition,  recent  analysis  based  on  a 
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limited  set  of  data  from  an  array  of  ionospherie  sounders  indicated  that  the  maximum 
density  of  the  anomaly  crest  and  the  latitudinal  width  of  the  anomaly  have  linear 
dependence  on  the  peak  value  of  the  vertical  plasma  drift  velocity  [Whalen,  2001], 

3.  Technical  Approach 

The  technical  approach  for  phase  1  is  to  calculate  the  electron  density 
distributions  as  a  function  of  altitude  and  latitude  at  a  specific  longitude  and  time,  using 
the  AFRL  theoretical  model  of  the  low-latitude  ambient  ionosphere.  The  model  uses  the 
following  input  parameters  as  a  function  of  location,  time,  solar  activity,  and 
geomagnetic  activity: 

1 .  The  neutral  densities  and  the  neutral  temperature  are  obtained  from  the  MSIS 
model  of  Hedin  [1986]. 

2.  The  horizontal  components  of  the  neutral  wind  velocity  are  obtained  from  the 
empirical  wind  model  HWM  of  Hedin  etal.  [1991]. 

3.  The  electron  and  ion  temperatures  are  based  on  the  empirical  model  of  Brace 
and  Theis  [1981], 

4.  Production  and  loss  rates  are  calculated  from  the  local  model  of  Jasperse 
[1974];  production  rates  use  the  solar  EUV  fluxes  derived  from  the  algorithm 
of  Hinteregger  etal.  [1981]. 

5.  Zonal  ion  drifts  are  given  by  an  empirical  model  [Fejer  et  al,  1991]  based  on 
Jicamarca  radar  measurements.  When  vertical  drift  velocity  information  is  not 
available,  the  model  uses  the  empirical  drift  model  [Scherliess  and  Fejer, 
1999]  based  on  measurements  by  the  Jicamarca  radar  and  the  AE-E  satellite. 
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Obviously,  if  real-time  data  for  neutral  densities,  neutral  wind  velocity,  electron 
and  ion  temperatures,  and  plasma  drift  velocities  were  available,  the  physics-based  model 
would  reproduce  the  low-latitude  ambient  ionosphere  quite  accurately.  Unfortimately, 
that  is  not  the  case  and  we  have  to  rely  on  models  for  these  important  input  parameters. 
The  adopted  models  for  the  input  parameters  do  a  fairly  good  job  in  reproducing  the 
climatology;  but,  for  the  day-to-day  variability  of  the  parameters,  improved  models  are 
desirable.  The  in-situ  and  ground-based  measurements  of  neutral  densities  and  neutral 
wind  velocity,  in  coordination  with  the  modeling  efforts  supplemental  to  the 
measurements,  such  as  TIGCM,  can  provide  the  desired  improvements.  The  C/NOFS 
mission  is  designed  for  such  a  purpose.  In  another  area  needing  improvement,  it  is 
necessary  to  have  either  measurements  or  a  physics-based  model,  instead  of  the  current 
empirical  model,  for  electron  and  ion  temperatures,  which  determine  both  the 
recombination  and  the  diffusion  rates.  The  empirical  model  of  Scherliess  and  Fejer 
{1999]  for  the  vertical  plasma  drift  velocity  does  a  fairly  good  job  in  reproducing  the 
climatology  of  the  temporal  behavior  of  the  plasma  drift  velocity.  But  it  often  fails  to 
predict  the  peak  value  of  the  post-sunset  velocity  accurately.  A  physics-based  model  is 
being  developed  at  AFRL  to  calculate  the  enhanced  eastward  electric  field  after  sunset.  In 
the  absence  of  drift  data,  this  model,  when  combined  with  the  empirical  model,  should 
yield  a  much  better  drift  model.  Alternatively,  if  the  peak  value  of  the  drift  velocity  can 
be  estimated  from  other  measurements,  as  the  present  study  suggests,  that  information 
may  also  be  incorporated  into  the  empirical  model  to  make  it  more  realistic. 

In  the  present  study,  vertical  plasma  drift  velocities  measured  by  the  Jicamarca 
incoherent-scatter  radar  (11.95°  S,  283.13°  E  geographic,  magnetic  dip  2°  N)  are  used  in 
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the  plasma  density  model.  Specifically,  F  region  average  values  of  the  drift  velocity  for 
fourteen  days  in  2002  (April  15-17,  May  31  -  June  1,  June  3-4,  October  8-10,  and 
November  11  -  14),  representing  both  equinox  and  solstice  conditions,  are  used.  These 
are  the  only  days  in  2002  for  which  the  Jicamarca  drift  data  are  available.  However,  the 
peak  values  of  the  post-sunset  drift  velocity  on  these  days  range  from  about  10  m/s  to 
about  70  m/s  and  so  the  data  set  is  useful  for  the  present  study.  Solar  conditions  were 
moderate  for  those  days  with  the  daily  value  of  Fio.?  cm  flux  between  170  and  200  units, 
and  the  3  hr  average  of  daily  Kp  values  ranges  fi'om  1.5  to  6.0. 

One  of  the  many  outputs  of  the  model  is  the  latitudinal  profile  of  the  height- 
integrated  electron  density  (TEC)  at  a  specified  longitude  (longitude  of  Jicamarca  radar  in 
the  present  study)  and  at  different  times.  For  validation  purpose,  the  latitudinal  profiles 
given  by  the  model  are  compared  with  those  obtained  fi'om  the  GPS  data  and  the  dual 
fi-equency  altimeter  data  fi-om  TOPography  Experiment  (TOPEX),  when  available.  The 
GPS  data  for  TEC  are  obtained  from  a  low-latitude  chain  of  10  receivers  [Valladares  et 
al,  2001]  that  are  located  near  the  west  coast  of  South  America  (-287°  E  longitude), 
which  is  close  to  the  longitude  of  the  Jicamarca  radar,  extending  between  5°  N  and  37°  S 
(geographic  latitudes).  From  the  latitudinal  distribution  of  the  calculated  TEC,  the  various 
anomaly  parameters  such  as  the  location  and  magnitude  of  the  anomaly  crest,  the  width 
of  the  anomaly,  and  the  peak-to-valley  ratio  of  anomaly  TEC,  etc.  are  determined. 
Repeating  the  calculations  of  TEC  for  the  fourteen  days  in  2002,  coefficients  of 
correlation  between  the  anomaly  parameters  and  the  peak  value  of  the  drift  velocity  are 
calculated  and  then  theoretical  relationships  between  them  are  determined  by  means  of 
least  squares  fits. 
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4.  Results 


The  four  panels,  (a)  -  (d),  in  Figure  1  show  examples  of  the  typical  outputs  of  the 
AFRL  ambient  ionosphere  model.  Figure  la  shows  the  plasma  drift  velocity  measured  by 
the  Jicamarca  radar  on  15  April  2002  as  a  function  of  universal  time  (UT)  [~  5  hours 
ahead  of  local  time  (LT)],  which  is  used  for  obtaining  the  results  presented  in  Figures  lb 
-  Id.  It  shows  the  daytime  enhancement  of  the  drift  velocity  near  -0700  UT  (~  1200  LT) 
due  to  the  E-region  dynamo  electric  field  and  the  post-sunset  stronger  enhancement  near 
0000  UT  (~  1900  LT)  due  to  the  F-region  dynamo  electric  field,  after  which  the  drift 
velocity  decreases  and  reverses  its  direction  around  0100  UT  (~  2000  LT).  Almost  the 
same  pattern  is  found  for  the  drift  velocity  on  all  fourteen  days.  Figure  lb  shows  the 
calculated  TEC  at  the  magnetic  latitude  (ML AT)  16“N  and  at  the  Jicamarca  longitude  as 
a  function  of  local  time.  The  magnetic  latitude  16°  N  has  been  chosen  because  the 
northern  crests  of  both  the  daytime  and  the  nighttime  anomaly  occur  near  it.  The 
maximum  TEC,  near  1600  LT  and  2100  LT,  represents  the  daytime  and  nighttime 
anomalies,  respectively.  Figure  Ic  shows  the  latitudinal  profiles  of  TEC  at  different 
selected  times.  It  shows  a  weak  (small  peak-to-valley  ratio  of  TEC)  anomaly  near  sunset 
(1800  LT).  After  sunset,  as  the  evening  progresses,  a  well-defined  anomaly  begins  to 
develop,  and  the  TEC  at  the  anomaly  crests  attains  its  maximum  value  around  2100  LT 
(shown  in  solid  line)  with  the  locations  of  the  anomaly  crests  at  about  ±16°  magnetic 
latitude.  At  later  times,  the  TEC  at  the  anomaly  crests  decreases,  the  locations  of  the 
anomaly  crests  move  toward  the  magnetic  equator,  and  eventually  in  the  morning  hours 
the  anomaly  disappears.  This  temporal  behavior  of  the  anomaly  is  in  good  agreement 
with  that  observed  experimentally.  Finally,  Figure  Id  shows  the  color-coded  map  of  TEC 
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Figure  1.  Example  of  the  outputs  of  the  AFRL  ambient  ionosphere  model,  (a)  Drift 
velocity  measured  by  Jicamarca  radar  on  15  April  2002  as  a  function  of  universal  time 
(UT).  (b)  Model  calculated  TEC  at  the  magnetic  latitude  (MEAT)  16°  N  and  at  the 
Jicamarca  longitude  as  a  function  of  local  time,  (c)  Latitudinal  profiles  of  calculated  TEC 
at  different  selected  times,  (d)  Color-coded  map  of  TEC  as  a  function  of  geographic 
latitude  and  universal  time. 
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as  a  function  of  geographic  latitude  and  universal  time  obtained  for  the  drift  data  of  April 
15.  It  clearly  shows  the  two  anomaly  crests  and  the  temporal  variation  of  their 
magnitudes,  as  well  as  the  TEC  minimum  near  the  magnetic  equator. 

Figures  2-5  show  comparisons  between  the  latitudinal  profiles  of  calculated 
TEC  and  those  obtained  from  GPS  measurements.  The  chosen  examples  represent 
different  months  in  2002.  The  four  panels,  (a)  -  (d),  in  each  figure  are  for  four  different 
times  indicated  therein.  The  solid  curve  in  each  panel  represents  GPS  data  and  the  dashed 
curve  represents  model  results.  The  dots  in  panel  (c)  of  Figure  3  and  panel  (b)  of  Figure  4 
represent  TOPEX  data.  Considering  the  fact  that  real-time  data  has  been  used  for  just  one 
of  the  several  input  parameters  of  the  model,  the  agreements  between  model  results  and 
observations  can  be  considered  to  be  reasonably  good. 

The  correlation  between  the  peak  value  of  the  post-sunset  drift  velocity  and  the 
various  anomaly  parameters,  such  as  the  location  and  magnitude  of  the  anomaly  crest,  the 
latitudinal  width  of  the  anomaly,  and  the  peak-to-valley  ratio  of  TEC  at  the  anomaly, 
obtained  from  the  theoretical  model,  have  been  studied.  It  is  found  that  significant 
correlation  exists  only  between  the  peak  value  of  the  post-sunset  drift  velocity  and  the 
peak-to-valley  ratio  of  anomaly  TEC.  The  six  panels,  (a)  -  (f),  in  Figure  6  show  the 
scatter  plots  of  the  peak-to-valley  ratio  of  model  calculated  TEC  in  the  northern  anomaly 
vs.  maximum  drift  velocity  (VM)  for  six  selected  times  indicated  therein.  The  straight 
line  in  each  panel  is  the  least  squares  fit  to  the  scatter  plot.  The  slope  and  intercept  of  the 
straight  line  as  well  as  the  value  of  the  correlation  coefficient  (R)  are  indicated  in  each 
panel.  At  0000  UT  (~  1900  LT),  which  is  about  the  time  when  vertical  plasma  velocity  on 
each  day  reaches  its  maximum  value,  the  peak-to-valley  ratios  are  almost  independent  of 
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Figure  2.  TEC  on  April  16,  2002  vs.  magnetic  latitude  at  four  selected  times:  (a)  0000 
UT;  (b)  0100  UT;  (c)  0200  UT;  (d)  0300  UT.  The  dashed  curves  represent  model  results 
and  the  solid  curves  represent  GPS  data. 
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Figure  3.  Same  as  in  Figure  2,  but  for  June  2,  2002.  The  dots  in  panel  (c)  represent 
TOPEX  data. 
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Figure  6.  Scatter  plots  of  the  peak-to-valley  ratio  of  model  calculated  TEC  in  the  northern 
anomaly  vs.  maximum  drift  velocity  (VM)  for  six  selected  times:  (a)  0000  UT;  (b)  0100 
UT;  (c)  0200  UT;  (d)  0300  UT;  (e)  0400  UT;  (f)  0500  UT.  The  straight  line  in  each  panel 
is  the  least  squares  fit  to  the  scatter  plot.  The  slope  and  intercept  of  the  straight  line  as 
well  as  the  value  of  the  correlation  coefficient  (R)  are  indicated  in  each  panel. 
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the  drift  velocity.  However,  a  definite  linear  dependence  with  a  significant  correlation 
coefficient  becomes  evident  starting  at  0100  UT  (~  2000  LT).  This  is  physically 
explainable  since  it  takes  at  least  an  hour  for  the  ionosphere  to  ‘feel’  the  effect  of  the 
post-sunset  enhanced  drift  velocity.  The  linear  dependence  persists  for  about  three  hours. 
However,  both  the  slope  and  the  intercept  of  the  regression  line  change  with  time.  This  is 
not  surprising  since  the  only  constant  feature  of  the  Earth’s  ionosphere  is  its  variability. 
The  six  panels,  (a)  -  (f),  in  Figure  7  show  the  scatter  plots  of  the  peak-to-valley  ratio  of 
the  model  calculated  TEC  in  the  southern  anomaly  vs.  maximum  drift  velocity  (VM)  for 
six  selected  times  indicated  therein.  As  in  the  case  of  the  northern  anomaly,  a  linear 
dependence  is  found  starting  at  0100  UT  (~  2000  LT)  and  seems  to  persist  for  a  couple  of 
hours.  However,  the  correlation  coefficients  for  the  southern  anomaly  are  evidently 
smaller  than  those  found  for  the  northern  anomaly.  The  six  panels,  (a)  -  (f),  in  Figure  8 
show  the  scatter  plots  of  the  location  of  the  northern  anomaly  crest  (LAC)  in  magnetic 
latitudes  vs.  the  maximum  drift  velocity  (VM)  for  the  six  selected  times  indicated  therein. 
The  straight  line  in  each  panel  is  the  least  squares  fit  to  the  scatter  plot.  The  slope  and 
intercept  of  the  straight  line  as  well  as  the  value  of  the  correlation  coefficient  (R)  are 
indicated  in  each  panel.  The  six  panels  (a)  -  (f),  in  Figure  9  show  similar  plots  for  the 
southern  anomaly  crest.  Note  that  the  correlation  coeffieients  for  the  location  of  the 
anomaly  erest  are  smaller  than  those  for  the  peak-to-valley  ratio  of  TEC.  The  coefficients 
of  correlation  for  the  other  anomaly  parameters  are  found  to  be  even  smaller.  In  Figures  6 
-  9,  results  for  the  days  when  the  peak  value  of  the  drift  veloeity  is  smaller  than  30  m/s 
have  been  excluded.  They  are  excluded  because  the  empirical  studies  of  Fejer  et  al. 
[1999]  and  Whalen  [2001]  indicate  a  threshold  value  of  the  peak  vertical  plasma  velocity 
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Figure  8.  Scatter  plots  of  the  location  of  the  northern  anomaly  crest  (LAC)  in  magnetic 
latitudes  vs.  the  maximum  drift  velocity  (VM)  for  six  selected  times:  (a)  0000  UT;  (b) 
0100  UT;  (c)  0200  UT;  (d)  0300  UT;  (e)  0400  UT;  (f)  0500  UT.  The  straight  line  in  each 
panel  is  the  least  squares  fit  to  the  scatter  plot.  The  slope  and  intercept  of  the  straight  line 
as  well  as  the  value  of  the  correlation  coefficient  (R)  are  indicated  in  each  panel. 
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for  the  scintillation  to  occur.  This  threshold  value  is  40  m/s. 


The  six  panels,  (a)  -  (f),  in  Figure  10  show  the  scatter  plots  of  the  peak-to-valley 
ratio  of  GPS  measured  TEC  in  the  northern  anomaly  vs.  maximum  drift  velocity  (VM) 
for  six  selected  times  indicated  therein.  The  straight  line  in  each  panel  is  the  least  squares 
fit  to  the  scatter  plot.  The  slope  and  intercept  of  the  straight  line  as  well  as  the  value  of 
the  correlation  coefficient  (R)  are  indicated  in  each  panel.  The  correlation  coefficients  are 
small  at  all  times.  The  six  panels,  (a)  -  (f),  in  Figure  1 1  show  the  peak-to-valley  ratios  of 
GPS  measured  TEC  in  the  southern  anomaly  vs.  maximum  drift  velocity  (VM)  for  six 
selected  times  indicated  therein.  The  correlation  coefficients  (R)  are  small  in  this  case 
too.  If  the  somewhat  stray  data  point  corresponding  to  VM  =  51  m/s  is  excluded  from 
Figure  11,  the  correlation  coefficients  at  0100  UT  and  0200  UT  become  comparable  to 
those  obtained  from  the  model  results  (see  Fig.  7).  This  is  shown  in  Figure  12.  Similar 
improvements  in  the  correlation  coefficients  are  not  obtained  if  we  exclude  the  data  point 
corresponding  to  VM  =  51  m/s  from  Figure  10. 

5.  Discussion 

Theoretical  study  of  the  nighttime  equatorial  anomaly,  based  on  the  low-latitude 
ambient  plasma  density  model  at  AFRL,  shows  that  linear  relationships,  with  signifieant 
correlation  coefficients,  exist  between  the  peak  value  of  the  post-sunset  enhanced  plasma 
drift  velocity  and  the  peak-to-valley  ratio  of  TEC  in  both  the  northern  and  the  southern 
anomaly.  The  relationships  are  valid  for  the  longitudinal  sector  of  the  Jicamarca 
incoherent  scatter  radar  whose  drift  velocity  measurements  have  been  used  in  the  model. 
Comparisons  of  model  results  with  available  data  show  reasonably  good  agreements  on 
certain  days,  but  not  on  other  days.  The  theoretical  calculations  use  real-time  data  for 
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Figure  10.  Scatter  plots  of  the  peak-to- valley  ratio  of  GPS  measured  TEC  in  the  northern 
anomaly  vs.  maximum  drift  velocity  (VM)  for  six  selected  times:  (a)  0000  UT;  (b)  0100 
UT;  (c)  0200  UT;  (d)  0300  UT;  (e)  0400  UT;  (f)  0500  UT.  The  straight  line  in  each  panel 
is  the  least  squares  fit  to  the  scatter  plot.  The  slope  and  intercept  of  the  straight  line  as 
well  as  the  value  of  the  correlation  coefficient  (R)  are  indicated  in  each  panel. 
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Figure  12.  Same  as  in  Figure  1 1  except  that  the  data  point  corresponding  to  VM  s  51  m/s 
has  been  excluded.  It  shows  improvement  in  the  values  of  R  at  0100  and  0200  UT. 
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only  one  input  parameter,  namely,  the  vertical  plasma  drift  velocity.  Empirical  models 
specify  other  input  parameters.  Either  real-time  data  or  improved  empirical  models,  such 
as  those  discussed  earlier,  for  the  input  parameters  need  to  be  used  in  order  to  have  more 
confidence  in  the  theoretical  model  results. 

The  theoretical  relations  presented  here  may  be  used  to  estimate  the  peak  value  of 
the  drift  velocity  from  the  measured  value  of  the  peak-to-valley  ratio  of  anomaly  TEC. 
The  time  dependence  of  the  relations  suggests  that  the  relation  that  is  appropriate  at  the 
time  of  measurement  should  be  used.  Such  measurements  are  provided  by  GUVI  on 
TIMED  satellite,  and  will  be  provided  by  UV  instruments  on  NPOESS  satellite.  In  the 
absence  of  drift  velocity  data,  the  estimated  peak  value  of  the  drift  velocity  may  be 
incorporated  into  the  empirical  drift  model  of  Scherliess  and  Fejer  [1999]  to  make  it 
more  realistic.  This  drift  model  can  then  be  used  in  the  AFRL  model  for  equatorial 
plasma  bubble  formation  to  forecast  the  occurrence  of  scintillation  and  the  magnitude  of 
the  amplitude  scintillation  index  S4. 

Two  important  issues  have  not  been  addressed.  One  important  issue  is  the  initial 
assessment  of  the  improvement  in  the  forecast  of  scintillation  offered  by  this  approach. 
The  other  important  issue  is  whether  or  not  similar  relationships  between  the  peak  value 
of  the  drift  velocity  and  the  anomaly  parameters  exist  in  other  longitudinal  sectors.  These 
issues  should  be  addressed  in  phase  2. 

After  the  work  was  completed,  it  came  to  oiu  attention  that  Bramley  and  Young 
[1968]  also  investigated  the  effects  of  drift  velocity  on  the  formation  and  the  structure  of 
the  equatorial  Appleton  anomaly  and  concluded  that  a  linear  relationship  exists  between 
the  drift  velocity  and  the  peak-to-valley  ratio  of  anomaly  TEC.  However,  their 
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calculations  and  analyses  were  done  for  steady-state  electron  density  distributions.  The 
steady-state  assumption  is  not  a  realistic  assumption  since  the  formation  of  the  anomaly 
is  a  dynamic  process  with  plasma  loss  and  diffusion  as  well  as  the  plasma  drift  velocity 
being  time-dependent.  The  results  of  our  calculations  done  with  the  time-dependent 
AFRL  ionosphere  model,  assimilating  drift  velocity  data,  should  be  much  more  accurate. 
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